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Abstract 



Backward proton and 7r~ production has been studied in CC interactions with 
carbon nuclei. Detailed analyses of the momentum distributions, of the production 
rates, and of the general features of events with a backward going particle, have been 
carried out in order to understand the mechanism producing these particles. The 
backward proton data have been compared with the predictions of the reinteraction 
and the short range correlation models. 

Key words: Neutrino interactions, cumulative production, intranuclear cascade, 
short range correlations 
PACS: 13.15.+g, 13.85.Ni 



1 Introduction 



It is a well established experimental fact that in the high energy interactions 
off nuclei there are particles emitted backwards, with respect to the beam di- 
rection, which have energies not allowed by the kinematics of collisions on a 
free and stationary nucleon. Backward going protons are commonly observed 
while, in absence of nuclear effects, their production is forbidden. High energy 
mesons, whose production in the backward direction is only allowed up to a 
given momentum, are detected also, at momenta above such a limit. Since a 
long time [1] this effect has been used as a powerful tool to investigate nu- 
clear structure. The models proposed to explain the origin of these particles 
(also called cumulative in the literature) can be divided essentially into two 
categories: models based on the intranuclear cascade mechanism and models 
based on the cumulative effect of groups of correlated nucleons/ quarks. 
In the intranuclear cascade models, the production of particles in the kinemat- 
ically forbidden region (KFR) can be seen as the result of multiple scattering 
and of interactions of secondary hadrons, produced in the primary zz-nucleon 
collision, with the other nucleons while they propagate through the nucleus 
[2]. The reinteraction or intranuclear cascade (INC) models usually rely on 
Monte Carlo methods to make their predictions [3], [4]. The importance of 
the intranuclear cascade mechanism is that it can provide information on the 
space-time evolution of the hadronization process. Experimentally one ob- 
serves that the cascade is restricted to slow particles only, while the fast ones 
do not reinteract inside the nucleus. The currently accepted explanation for 
this effect is the "formation zone" concept [5] -[7]. This is the distance (or the 

1 Now at Univ. of Perugia and INFN, Perugia, Italy 

2 Corresponding author. Email address: veltri@fis.uniurb.it 
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time) from the production point which is required for the secondary hadrons 
to be "formed", i.e. to be able to interact as physical hadronic states. Since 
the distance/time required, due to the Lorentz time-dilation factor, is pro- 
portional to the energy of the secondary, the INC process is restricted to slow 
hadrons which have formation lengths smaller than the nuclear radius. 
An advantage of neutrino (or more generally lepton) induced interactions with 
respect to hadronic processes is the fact that the projectile interacts only once, 
avoiding the complications related to the projectile rescattering in the target. 
Natural drawbacks are the facts that the equivalent of the projectile energy in 
the hadronic processes, i.e. the hadronic jet energy in the leptonic processes, 
is not fixed but varies from event to event and suffers from systematic uncer- 
tainties, especially related to the reconstruction of neutrals. 
Another feature of the neutrino-nucleus interaction is the fact that, due to 
the extremely small neutrino-nucleon cross-section, the interaction can practi- 
cally take place anywhere in the nucleus. This is not the case in hadron-nucleus 
(and also photon-nucleus) interactions where, due to the large cross-section, 
the interaction takes place essentially on the first nucleons along the particle 
path, i.e. on the surface of the nucleus. As a consequence not all the nucle- 
ons participate in the scattering and the total hadron-nucleus cross-section 
is smaller than A (the number of nucleons in the nucleus) times the total 
hadron-nucleon cross-section. This effect is known as shadowing [8]. 
In the correlated nucleon/ quark models the backward particles are produced 
in the collisions off structures with mass larger than the mass of the nucleon. 
These structures are formed, at small interparticle distance, under the action 
of the short range part of the nuclear force. They may either be described as 
fluctuations of the nuclear density [9] or as clusters of a few correlated nucleons 
[10]- [12]. In any case these structures represent the effect of gathering two or 
more nucleons in small volumes with a radius of the order of 0.5^-0.7 fm. The 
nucleons in these structures can acquire high momenta and the fast backward 
going particles can be seen as a direct manifestation of the high momentum 
tail of the Fermi distribution. For instance, the classical mechanism of Fermi 
motion of nucleons inside the nucleus can explain backward proton produc- 
tion only up to about 300 MeV/c, well below the observed limit. Nucleons can 
also loose their identity in these structures and larger masses can be reached 
if quarks in neighbouring nucleons stick together, forming multiquark clusters 
[13,14]. The momentum distribution of quarks in the cluster is, in this case, 
responsible for the observed spectra of cumulative particles. 
These two classes of models are not mutually exclusive. Indeed some features 
of the data can be explained by both, and the experimentally observed pro- 
duction in the backward hemisphere [15] can have contributions from both 
mechanisms. 

While the production of particles in the KFR has been investigated over a very 
wide range of incident energies and on many different nuclear targets using 
hadron beams [16], the data on backward particle production in v (&) induced 
reactions are far less abundant. Up to now only five bubble chamber experi- 
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merits have studied backward protons (and in only one experiment backward 
pions as well) using bubble chambers filled with deuterium [20], heavy liquid 
fillings [17]- [20] and in one case a hybrid emulsion-bubble chamber technique 
[21]. " ^ 

Early hadronic experiments had found that the invariant cross-section for 
backward particle production can be parametrized as e~ BP where P is the 
particle momentum. It was also found that the slope parameter B is almost 
independent of the type of incident particle, its energy and target type, a fact 
known as "nuclear scaling" . Neutrino experiments have confirmed these prop- 
erties. 

Here we present a detailed study of backward going protons {Bp ) and back- 
ward going 7r~ (Bit~ ) produced in charged current neutrino interactions in 
the NOMAD detector. In Sec. 2 and Sec. 3 we describe the experimental ap- 
paratus and the data treatment. In Sec. 4 we report the invariant cross section 
distributions of Bp and Btc" . In Sec. 5 we discuss the kinematical properties 
of events with a backward particle. In Sec. 6 we give the production rates 
and study their dependence on the atomic number using the results of other 
neutrino and hadron experiments as well. Finally in Sec. 7 we compare the 
data with the predictions of theoretical models. 



2 NOMAD detector and neutrino beam 

For the study of backward particles presented here, the most important com- 
ponent of the NOMAD detector is the active target. Its main features are 
highlighted here, while a detailed description of the full detector (shown in 
Fig. 1) can be found in Ref. [22] and [23]. 

The target was designed with two conflicting requirements. It had to be as 
light as possible in order to minimize photon conversions, multiple scattering 
and secondary hadronic interactions, and it had to be as heavy as possible in 
order to produce a large number of neutrino interactions. 
This conflict was resolved using an active target (2.7 tons) of 44 drift chambers 
(DC), 3 x 3m 2 each [24], perpendicular to the beam axis. The target mass is 
provided by the chamber structure having an average density of 0.1 g/cm 3 . 
The drift chamber composition is reported in Table 1, where it can be seen 
that carbon and elements with nearby atomic numbers represent over 90% 
of the total weight. For this reason we consider this study as a measurement 
of backward production in v^C interactions. The chambers are placed inside 
a 0.4 T magnetic field and provide a momentum resolution which can be 
parametrized as: 

op _ 0.05 0.008P 
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where L is the track length in meters and P the momentum in GeV/c. The 
first term is the contribution from multiple scattering and the second comes 
from the single hit resolution. 
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Fig. 1. Side view of the NOMAD detector. 



Element 


Z 


Weight (%) 


p(%) 


n (%) 


H/D 


1 


5.14 


5.09 


0.05 


C 


6 


64.30 


32.12 


32.18 


N 


7 


5.92 


2.96 


2.96 





8 


22.13 


11.07 


11.07 


Al 


13 


1.71 


0.82 


0.89 


Si 


14 


0.27 


0.13 


0.14 


CI 


17 


0.30 


0.14 


0.16 


Ar 


18 


0.19 


0.09 


0.10 


Cu 


29 


0.03 


0.01 


0.02 








52.43 


47.56 



Table 1 

NOMAD drift chamber composition, showing the proportions by weight of atoms, 
and of protons and neutrons. As can be seen the target consists mainly of carbon 
and is practically isoscalar. 

The target is followed by a transition radiation detector (TRD) [25], a preshower 
(PRS) and an electromagnetic calorimeter (ECAL) [26]. A hadron calorimeter 
and two muon stations are located just after the magnet coil. The neutrino 
interaction trigger [27] consists of a coincidence between two planes of coun- 
ters located downstream of the active target, in the absence of a signal from 
a large area system of veto counters upstream of the NOMAD detector. 
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The CERN-SPS wide band neutrino beam is produced by 450 GeV/c protons 
incident on a beryllium target. Neutrinos are produced in the decay of sec- 
ondary pions and kaons in a 290 m long decay tunnel at an average distance 
of 625 m from the detector. The relative beam composition is predicted to 
be : : v e : v e = 1.00 : 0.0612 : 0.0094 : 0.0024 with average energies of 
23.5, 19.2, 37.1 and 31.3 GeV respectively [28]. The average energy of muon 
neutrinos interacting in the apparatus is about 41 GeV. 
Referring to Fig. 1 the coordinate system adopted for NOMAD has the X-axis 
into the plane of the figure, the F-axis upwards and the Z-axis horizontal, ap- 
proximately along the direction of the neutrino beam, which points upward, 
at an angle of 2.4° with respect to the Z-axis. Its origin is in the center of the 
front face of the first DC along the beam. 



3 The data sample 

3.1 Event selection 

This study is based on the NOMAD full data sample collected between 1995 
and 1998. In this analysis the event selection requires a primary vertex with 
at least 2 tracks inside the fiducial area defined by \X\ < 130 cm and — 125 < 
Y < 135 cm. Along the Z direction a cut Zvmin < Z < 400 cm is imposed. 
Zvmin is 5 cm for the 1995 and 1996 data and is 35 cm for the data collected 
in 1997 and 1998. In the second period the first DC module was removed to 
install NOMAD-STAR [29] and the fiducial volume was slightly reduced. 
To select a CC event a negative \x attached to the primary vertex, identi- 
fied by the muon stations and having a momentum of at least 3 GeV j c, was 
required. Under these conditions the total CC sample consists of 944019 
events. 

3.2 Track selection 

Only tracks attached to the primary vertex are used in the search for backward 
particles. The track is required to have at least 8 DC hits, the distance of the 
first hit from the primary vertex to be less than 15 cm, and the relative error 
on the reconstructed momentum to be less than 0.3. For the Bp search, since 
the identification method makes use of the momentum-range relation, tracks 
are required to range out in the DC volume. To avoid escaping tracks we 
require in addition the last hit of the track to be in a reduced fiducial volume 
defined by \X\ < 100 cm, \Y\ < 100 cm and Zmin < Z < 375 cm, where 
Zmin = 30 cm for the 1995 and 1996 data sets and Z MIN = 60 cm for the 
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rest. Finally no secondary vertex must be found close to the last hit of the 
track. 

3.3 MC sample 

In order to estimate proton and pion reconstruction efficiencies and back- 
ground we use a large sample of CC Monte Carlo (MC) events. These 
events are generated using a modified version of LEPTO 6.1 [30] and JETSET 

7.4 [31]. The target nucleon motion is simulated using the Fermi momentum 
distribution proposed by Bodek and Ritchie [33]. The generated events un- 
dergo a full detector simulation based on GEANT [32] and are subsequently 
reconstructed. MC events and tracks are selected in the same way as are ex- 
perimental data yielding a final MC sample consisting of ~ 2.5 x 10 6 events. 
The MC data used for this analysis do not include the simulation of "nuclear 
effects" such as the intranuclear cascade or correlations [23]. When needed, 
specific correction procedures to the MC events are applied using the experi- 
mental data distributions (see Sec. 3.5). 

3.4 Backward proton identification 

Fig. 2 shows the experimental distributions of length vs. momentum for pos- 
itive (left) and negative (right) tracks, going backward with respect to the 
beam direction, and satisfying our selection criteria. 




0.2 0.4 0.6 0.8 0.2 0.4 0.6 



P REC (GeV/c) P REC (GeV/c) 

Fig. 2. Distributions of length vs. momentum for positive (left) and negative (right) 
backward going tracks selected as described in the text. The line indicates the position 
of the cut (Eq. 1). 

Two distinct populations are clearly visible on the positive sample plot. Pro- 
tons, having a shorter range than 7r + , tend to accumulate in the lower right 
part of the plot while the 7r + 's tend to populate the left-hand side. Comparing 
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the two plots we see that the lower right part of the negative sample is much 
less populated than the corresponding region of the positive one, these tracks 
being mainly tt~. The separation between protons and pions is also visible in 
Fig. 3 where the momentum distribution of positive (dots) and negative (solid 
line) backward tracks is shown for different intervals of track length. In this 
figure the heights of the negative distributions have been normalized to the 
positive pion peaks. 

We identify as a proton any positive backward going track which passes our 
selection cuts and has length L: 

L < 2000 (P REC + 0.150) 3 6 cm (1) 



where Prec is the reconstructed track momentum in GeV/ c. The cut position 
was optimized by using the MC results of forward tracks and is shown in 
Fig. 2. 
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P BEC (GeV/c) 



Fig. 3. Reconstructed momentum distributions for positive (dots) and negative (solid 
line ) backward going tracks for different intervals of the track length. The heights of 
the negative distributions have been normalized to the positive pion peaks. 
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3. 5 Data corrections for the Bp analysis 

The reconstruction procedure used the energy loss of a pion for the track fit. 
Instead of refitting proton candidate tracks with the correct energy loss, an em- 
pirical correction was applied to the reconstructed proton momentum Prec- 
This correction, which was obtained from MC events, amounts to +180 MeV/c 
at Prec — 200 MeV/c and becomes negligible above 500 MeV/c. This pro- 
cedure gives the true proton momentum P. 

The raw number iV™ w of identified protons in a given momentum bin must be 
corrected for reconstruction efficiency ^rec-, stopping efficiency ^stop-, iden- 
tification efficiency e^o and pion contamination in the data ntont ta ■ All these 
corrections will be described below. They are functions of either the recon- 
structed or the true proton momentum. Since the MC used for this analysis 
does not properly account for Bp production, some of these efficiencies can 
only be obtained by correcting the MC efficiencies with the data, or directly 
from data themselves. The "true" number N p of protons is then: 

Np = N ra W J_ _J_ J_ (1 _ ^Jat a) (2) 

£rec (-stop £id 
with an average overall correction factor of 3.9 . 

a) Reconstruction efficiency £rec 

This is the ratio between the number of reconstructed and generated protons in 
a specific bin of the true momentum P and angle 9 (measured with respect to 
the beam direction). Since this quantity depends only on the detector geometry 
and composition, we assume that the predictions obtained using forward going 
protons are applicable to the backward ones with the replacement cos 9 — > 
— cos#. An average value of 0.48 is obtained for e RE c, the distribution of 
which is shown on the left in Fig. 4. 

b) Stopping efficiency €stop 

This is the ratio between the number of backward protons fully contained in 
the target (Bpsrop) and all the reconstructed ones (Bp ). Its average value is 
0.42 in the momentum range used in this analysis. 

The protons are extracted from the sample of positive backward tracks by 
subtracting the backward n + content. The stopping efficiency esrop is com- 
puted from the following quantities: 

BpsTOP _ B$top ~ Bpstop 
esroP - — ^- - B+ _ Bn+ 

where B + (Bg TOP ) refer to backward (backward contained) positive tracks. 
Bn + and Bng TOP are the analogous quantities for backward positive pions. 
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Fig. 4. Reconstruction efficiency (as described in the text) as a function of momen- 
tum and angle between the track and the nominal beam direction. The plot on the 
left is for protons, the one on the right for ir~ . 

They are evaluated from the data as follows. 

In the NOMAD isoscalar target the length distributions of ir + and n~, in each 
momentum bin, are the same apart from a scale factor. Therefore the number 
of backward 7r + can be computed by the number of backward 7r~ once the 
relative population of Bn + and Bti~ is known. The ratio of the two popula- 
tions is measured from a clean sample of pions obtained by selecting tracks 
with a sufficiently longer length than the corresponding proton range. This 
ratio has comparable values for contained and not contained pions. It is al- 
most constant in the momentum range of interest with a weighted average of 
R = 1.82 ±0.05. 

The pion content in the backward sample is therefore: 

Bn + = R x Bn~ ; Birg TOP = R x Biig TOP 

where the negative pion sample (Bti~ , B^g TOP ) is taken to be the sample of 
negative (negative contained) backward tracks. This relies on the fact that the 
MC simulation shows that the e ± contamination in the backward region is at 
the level of a few percent above 200 MeV j c. 

c) Identification efficiency ejo 

This is the fraction of Bp with track length smaller than the length cut. Above 
Prec — 300 MeV jc this quantity is 1. The results of forward going protons 
are applicable also in this case. 

d) Pion contamination T^tont ta 

This is the fraction of 7r + 's in the sample of identified protons, which amounts 
to ~ 8% above Prec — 250 MeV jc. To estimate this contamination in the 
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data, the MC prediction (7r^ ni ) for forward particles has been corrected us- 
ing the backward data distributions as follows. For each momentum bin the 
length distributions of protons in the data and MC can differ in normalization 
but not in shape since the track length depends on the mechanism of energy 
loss which is well reproduced by the MC. The same applies to 7r + . The length 
cut applied on the momentum-length plane (Eq. 1) defines in the data and 
MC two populations: the 7r-like tracks, (the tracks above the cut) and the 
p-like tracks (those below the cut). To account for the difference in the popu- 
lations of protons and pions in the data and MC the quantity ^tont° has been 
weighted by the double ratio of the 7r-like tracks (data/MC) and p-like tracks 
(data/MC): 

]\jdata Tvrdata I j\tMC 

+ data . . iV 7r+ . _ tt — like I n — like + MC 
n cont ~ pjdata ~ jydata /jyMC cont 

The above correction is evaluated as a function of the reconstructed momen- 
tum. 



3. 6 Backward tt identification 

As already pointed out the e~ component is negligible in the negative back- 
ward track sample. We therefore assume that any negatively, backward going, 
charged track with P > 0.2 GeV/c is a tt". The smallness of the contamination 
and the fact that it is not necessary to look for stopping tracks to identify the 
Btt~ reduce the number of corrections to be applied to the data to only one, 
the reconstruction efficiency e^j^.For this quantity we use the MC predictions 
obtained for forward going 7r~. 

The true number of negative pions is then obtained from the raw number of 
identified negative pions jV™ w in a given momentum bin according to: 

N*- = K aw 4- ( 3 ) 

e REC 

In this case the average value of the correction is 2. 



4 Backward p and it invariant momentum distributions 

The inclusive spectrum of backward particles is typically represented using 
the normalized invariant cross section (1/<7tot) E d 3 a/dP 3 , where E is the 
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energy of the backward going particle. The invariant cross section is usually 
[20] parametrized by an exponential form as: 



1 E dN 
NT v Pdp5 



= Ce 



-BP Z 



(4) 



where N ev is the total number of events, C is a constant and B is the slope 
parameter. In previous experiments B has been found to be almost indepen- 
dent of the projectile type and momentum, and of the atomic number of the 
target. This behaviour has been termed "nuclear scaling [16]. 
The inclusive spectrum of Bp and Btt~ is shown in Fig. 5a and 5b respectively, 
together with the exponential fit. The values of the measured slope parameter 
B are reported in Table 2. The first and the second errors are statistical and 
systematical, respectively. 
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Fig. 5. Invariant momentum distributions for backward going protons (a), ir (b) 
in v^CC events and protons in v^CC events (c) 



The systematic uncertainty was estimated by slightly changing the values of 
all the cuts used, by varying by a small amount the correction functions and 
also by changing the size of the fiducial volume. The values of the slope pa- 
rameters measured in this experiment are found to be compatible with the 
results obtained in other neutrino (see Fig. 6 for the Bp case) and hadronic 
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ev i c) 


ij^c /uev ) 


o (c /Lrev ) 


Backward p 




CC 


0.37^ 


0.70 


10.54 ±0.20 ±0.5 


4.08 ±0.19 ±0.5 






CC 


0.37 


0.70 


10.79 ± 0.78 


2.71 ±0.54 


Backward 7r~ 




CC 


0.32 


0.85 


10.03 ±0.28 ±0.3 


0.17 ±0.01 ±0.02 



Table 2 

Fit ranges and values of the slope parameter B, for backward protons and ir~ , as 
obtained from the exponential fit to the invariant cross section. The first and the 
second errors are statistical and systematical, respectively. 

experiments. The invariant cross section for Bp is larger than the one for 
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NOMAD (C) 



BEBC (D 2 ) 
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E-564 

(Em) 



Fig. 6. The value of the slope parameter B obtained in neutrino experiments. For 
NOMAD the statistical and systematical error have been added in quadrature. The 
full circles are the results from CC, the open circles from CC and the cross 
is for NC interactions. For BEBC v and 9 results have been combined. 

Bn~ by about one order of magnitude but the values of the slopes are similar. 
The kinematic ranges of the two invariant distributions are also different. To 
be identified Bp have to stop inside the target volume; the rather low density 
of the NOMAD target (see Sec. 2) restricts to m 0.7 GeV/c the maximum 
useful momentum value. 

The invariant cross section and slope parameter for Bp in CC events are 
given in Fig. 5c and Table 2 respectively. During normal operations CC events 
were also collected due to the small component of the dominantly beam 
[22]. A dedicated v run yielded an additional sample of CC events included 
for this analysis. 
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Antineutrino events are selected under the assumptions that the efficiencies 
and pion contamination are the same as those used for the neutrino events, 
but requiring a positively charged muon instead of a negative one. The final 
sample consists of 61134 events containing 1764 Bp . 

4-1 Energy dependence of the slope parameter 

In Fig. 7 we show the slope parameter B for protons as a function of the 
visible hadronic energy Eh ad and the Q 2 of the event. The visible hadronic 
energy is defined as: 

Eh ad = E c + E n 

where Y,E C is the sum of the energies of reconstructed charged tracks (as- 
suming the mass of the pion if the particle type is not explicitly identified). 
The sum includes all the charged tracks attached to the primary vertex by the 
reconstruction program, those belonging to a pointing V°, the "close hang- 
ers" i.e. isolated tracks not associated to the primary nor to a secondary 
vertex and having their first hit in a box around the primary vertex of size 
\X\, \Y\ < 10 cm and —5 < Z < 20 cm, and the "pointing hangers" i.e. 
the hangers which have an impact parameter of less than 10 cm when lin- 
early extrapolated back to the z-plane of the primary vertex. J2 E n includes 
tracks associated to secondary vertices corresponding to interactions of neu- 
tral particles and the energy of neutral particles reconstructed in the ECAL. 
The reconstructed neutrino energy E vis is taken to be the sum of the muon 
energy E^ and of Eh ad- The square of the four-momentum transfer Q 2 is 
Q 2 = 4 E vis E^sin 2 8/2, where 9 is the muon angle with respect to the neu- 
trino direction. No significant dependence of the slope B on either quantity is 
observed, in agreement with the expectations of "nuclear scaling" as observed 
in hadronic beam experiments. The range covered by NOMAD is similar to 
the one covered by different experiments with hadronic beams. 

4-2 Angular dependence of the slope parameter 

The angular dependence of the slope parameter B is usually given with respect 
to the direction of the exchanged boson, approximated by the hadronic jet 
direction [20]. This allows a comparison with the results obtained in hadron 
beam experiments. The angular dependence of B for Bp and Btc~ is shown in 
Fig. 8 as a function of the angles with respect to the beam and to the hadronic 
jet direction. Although these two angles are highly correlated, there are small 
differences in the slope values. If the beam direction is used B is systematically 
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Fig. 7. The value of the slope parameter B shown as a function of the hadronic 
energy Eh ad (top) and of Q 2 (bottom). 

smaller by rs 0.5^-l(Gel / /c)~ 2 , in the Bp case, while for B-ir~ the difference is 
less pronounced. As can be seen from Fig. 8, larger values of B are preferred at 
increasingly backward directions. This behaviour has already been observed 
in neutrino [20], photon [34] and hadron [35-37] nucleus experiments. 



5 Kinematical properties of events with a backward particle 

Table 3 shows a comparison of the average values of some kinematical quanti- 
ties for CC events with and without Bp or Btt~ . < C > and < N^ 1 > are 
the average event charge and charged multiplicity included), v = Emg—E^, 
x = Q 2 /2Mu is the Bjorken scaling variable and W 2 = M 2 — Q 2 + 2Mv is 
the square of the invariant mass of the hadronic system, M being the nucleon 
mass. 

For events with a backward particle the average values of the charged mul- 
tiplicity are consistent with the presence of extra tracks as expected in both 
models. The average value of the event charge in the Bp sample (< C >= 
2.11) is larger than in the no Bp, no Bn~ samples as measured in the data 
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Fig. 8. Angular dependence of the slope parameter B forp (left) and ir~ (right). The 
value of B is shown as a function of the particle angle with respect to the hadronic 
jet direction (full circles) and with respect to the beam direction (open circles). 





CC (no Bp, no Bit') 


CC (Bp) 


CC (Bir~) 


<c> 


0.46 


2.11±0.02 


0.08±0.02 


< N± > 


5.13 


6.63±0.03 


7.28±0.03 




40.75 


38.9±0.4 


39.5±0.5 


< Eh ad > 


12.32 


11.4±0.2 


11.7±0.2 


<Q 2 > 


7.07 


5.5±0.1 


5.9±0.1 


<x> 


0.32 


0.26±0.002 


0.28±0.003 


< W 2 > 


17.19 


16.4±0.2 


17.5±0.3 



Table 3 

Comparison of average values of some kinematical variables in events with and 
without an identified Bp or Bir~ . The values of < C > and of < > are 
corrected for the track reconstruction efficiency. Only statistical errors are shown. 
They are negligible for the first data column. 

(< C >= 0.46). It is also larger than in the total MC sample where no nu- 
clear effects are present (< C >= 0.31). This observed increase is not entirely 
due to the bias introduced by demanding a positive particle, since an average 
value of < C >= 1.21 is obtained when requiring a backward positive track 
in the MC events. This can be expected from both mechanisms outlined in 
the Introduction. In the collisions of secondary hadrons with nucleons inside 
the nucleus the total event charge will increase by one unit in interactions off 
protons while no extra charge will be produced when scattering on neutrons. 
Also in the framework of the short range correlation model (in its simplest 
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form of two-nucleon correlations) the mechanism of pair breaking predicts ex- 
tra tracks in the final state. The mechanism is of course symmetric for p and 
n but backward going neutrons are not detected. 

The events with a Btc~ show an average charge smaller by ~ 0.4 units with 
respect to the value where no Bn~ is present. Here again this effect cannot be 
entirely interpreted as the result of the bias associated with the requirement 
of at least one negative backward track, since the MC in this case predicts 
< C >= —0.41. Therefore also in the sample of events with a B7r~ an in- 
crease in the average value of the event charge is observed. 
From Table 3 it appears that the average values of the kinematical variables, 
other than < C > and < N^ 1 >, for events with either a Bp or a Bir~ are 
systematically lower than for events without identified backward particles. 



6 Backward particle rates 



The Bp rate has been compared with the results of the other //-nucleus ex- 
periments in order to study a possible atomic number dependence. In NO- 
MAD, where Bp are identified in the momentum interval from 370MeV/c 
to 700MeV/c we obtain, after correction, the number of events reported in 
Table 4. 



N. of Bp or Btt- 


N. of Bp ev. 


N. of Btt~ ev. 





904212 


939617 


1 


37634 


4238 


2 


2168 


164 


3 


5 






Table 4 

The number of events, after correction for efficiencies, as a function of the mul- 
tiplicity of backward protons in the momentum interval 370 700MeV/c, and of 
backward negative pions in the momentum interval 350 800MeV/c. 

These figures correspond to an average number of Bp per event, < Nb p >= 
(44.5 ± 0.5) x 10~ 3 , where the error is statistical only. The backward proton 
yield obtained in other //-nucleus experiments is shown in Table 5. The values, 
which were extracted from the original references, are not directly comparable 
because of the different momentum intervals used in the different experiments. 
In order to study the A dependence the Bp yields measured in NOMAD and 
E-564 [21] were extrapolated to the interval from 350 to 800 MeV/c assuming 
the dependence given in Eq. 4 with the measured slopes. 
The extrapolation for NOMAD gives: 



< N Rn > 



Bp ^350H-800MeV7c 



52.8 ± 0.6(stat.) ± 7(syst.) 



x 10 



-3 
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Experiment 


< A > 


A P 


<E V > 


< N Bp > 


< N Bn - > 






(GeV/c) 


(GeV) 


(xlO" 3 ) 


(xlO- 3 ) 


BEBC D 2 [38] 


2 


0.35^-0.80 


50 


6.3 ± 0.7 


1.0 ± 0.3 


NOMAD 


12 


0.37^-0.70 


41 


44.5 ± 0.5 


4.8 ± 0.1 


BEBC Ne [20] 


20 


0.35^0.80 


50 


97.7 ± 3.7 


6.8 ± 0.9 


15foot [17] 


20 


0.20^0.80 


50 


137.5 ± 15 




SKAT [18] 


52 


0.32^0.70 


10 


185.8 ± 10 




E-564 [21] 


80 


0.30^0.80 


60 


331.1 ± 47 





Table 5 

Comparison of Bp and Bn~ rates for neutrino-nucleus experiments. The average 
mass number < A >, the momentum interval AP used in the analysis, the average 
event energy < E u > and the Bp and Btt~ rates are listed for each experiment. The 
errors are purely statistical. 

For E-564 we obtain < N Bp >= (234 ± 33) x 1(T 3 . 

The A dependence for experiments most directly comparable to NOMAD is 
shown in Fig. 9a. In the range A = 20 4- 80 it has been parametrized as 
< N B p >oc A a , where a = 0.68 ± 0.12 [21]. The same parametrization with a 
similar power law was found to describe Bp production in n + and K + collisions 
with Al and Au nuclei at 250 GeV/c [39]. It is evident from Fig. 9a that this 
simple power law does not describe the NOMAD data taken at a lower < A >. 
The B7r~ rate was directly measured in the same momentum range used for 
the Bp analysis. Its average value is found to be: 



< N Bn - >350-800MeV/ c = 4.8 ± 0.1(stat.) ± 0.3(syst.) 



x 10 



-3 



In this case a good fit to the two BEBC points and to the NOMAD value can 
be obtained using the form A a giving a = 0.83 ± 0.25, as shown in Fig. 9b. 



7 Comparison of data with models 

7. 1 E H ad and Q 2 dependence of Bp and Bn~ rates 



The Bp and the Btt~ rates have been studied as a function of the hadronic 
energy E H ad and of Q 2 . In both cases, shown in Fig. 10 and 11 a decrease 
of the yield with increasing Ehad and Q 2 is observed. As pointed out in the 
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Fig. 9. The average number of Bp per event (a) and of Bn~ per event (b) in the 
momentum range from 350 800MeV/c in neutrino experiments as a function of 
the mass number A. The line shown in (b) is the result of the fit described in the 
text. 

Introduction, this can be interpreted in terms of the "formation zone" concept. 
The larger Eh ad and Q 2 , the larger is the average energy of the outgoing 
partons therefore resulting in hadrons which have higher probability to be 
formed outside the nucleus. As a consequence, reinteractions will decrease 
and so will the slow proton rates. This E HAD and Q 2 dependence is also 
consistent with the decrease of the average values of E HAD and Q 2 in events 
with identified Bp as shown in Table 3. In the B7i~ case the dependence of the 
yield on Q 2 and on Ehad is less pronounced. This can be understood since the 
Bit~ rates are a less sensitive probe of nuclear effects because Bn~ production 
on a stationary nucleon is kinematically allowed for momenta up to about half 
the nucleon mass. Furthermore Bn~ can be produced in the decay of forward 
going resonances or unstable particles. 

7.2 Multiplicities of slow particles 

Since mostly slow hadrons are produced inside the nucleus, they are predomi- 
nantly the ones that can scatter off nucleons, producing both slow protons and 
backward particles. Therefore, in the framework of the rescattering model, a 
correlation between the number of emitted slow protons and the multiplicity 
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Fig. 10. TTie average number of Bp (370 < P < 700MeV/c) per event as a function 
of the hadronic energy Eh AD (left) and of Q 2 (right). 
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Fig. 11. The average number of Bir~ (350 < P < SOOMeV/c) per event as a function 
of the hadronic energy Eh ad (left) a-nd of Q 2 (right). 

of mesons is expected [40]. 

The increase in the rate of events with slow protons as a function of hadron 
multiplicities was first observed by the E745 collaboration [41] in a z/-Freon 
experiment. This correlation was used [40] to extract information on the for- 
mation zone of secondary hadrons. This effect was also observed by the BEBC 
collaboration [42] . 

We have observed this correlation by studying the fraction of events with 
at least one Bp as a function of the multiplicity of low momentum (P < 
700MeV/c) positively or negatively charged hadrons (Fig. 12). In both cases, 
this fraction of events increases with increasing multiplicity of low momentum 
tracks. However, the positive hadronic multiplicity is biased by the presence 
of protons from rescattering itself and therefore it is less representative of 
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the original multiplicity of positive slow hadrons from the neutrino interac- 
tion. The true correlation can therefore be better studied as a function of the 
multiplicity of negative slow hadrons. 
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Fig. 12. Fraction Fb p of events with at least one Bp as a function of the multiplicity 
of low-momentum (P < 7Q0MeV/c) positively (left) or negatively (right) charged 
hadrons. 



7. 3 Backward tt spectra and the Fermi momentum tail 



As already pointed out in the Introduction, in the short range correlation 
models the spectrum of fast backward going hadrons reflects the tail of the 
Fermi momentum distribution. For this study we used Bir~ because, as op- 
posed to Bp, they do not need to range out in the target thus yielding data 
up to larger momenta. 

To estimate qualitatively the contribution of this tail we have compared our 
Bn~ data to the predictions of two different Fermi momentum distributions 
[33,43]. The distribution of Ref. [33] has a long momentum tail up to 4 GeV/c. 
In Ref. [43] the spectrum of backward going hadrons in lepton-nucleon scat- 
tering is derived from the spectral function p(k, E) , which is the probability 
of finding a nucleon with momentum k and removal energy E. The function, 
whose integral over E is the Fermi distribution, is obtained from non rela- 
tivistic many body theory and a specific procedure has been developed to 
extrapolate it to large values of k. In this case the Fermi momentum tail is 
shorter and ends at about 1 GeV/ c [44] . 

In Fig. 13 the invariant spectrum of MC Bir~, simulated using the Fermi dis- 
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P 2 (GeV 2 /c 2 ) 

Fig. 13. Invariant spectrum for Bir~ in MC (open circles) and data (full circles). 
The solid line is the curve of Ref. [43] superimposed with an arbitrary normalization. 

tribution of Ref. [33], is shown together with the data. A clear disagreement 
is visible in the tail of the distributions, the MC being much larger than the 
data above ~ 0.5 (GeV/c) 2 . Our results therefore do not support a Fermi mo- 
mentum distribution with a long tail as proposed in Ref. [33] while they agree 
with the dependence predicted by Ref. [43], as also shown in Fig. 13. 



7.4 Effects of short range correlations in Bp production 



According to the picture proposed in Ref. [10] Bp production is explained as 
the result of a neutrino interaction within a correlated cluster of two or more 
nucleons. This cluster is formed, for a short time, when two nucleons in their 
motion inside the nucleus approach each other so closely as to come under the 
effect of the short range component of the nuclear force (r c = 0.5 4- 0.7 fm). 
As a result, the high relative momenta of the correlated nucleon pair manifest 
themselves when the backward moving spectator is released in the interaction 
of the incoming virtual W with the forward going nucleon. In this model, if 
the effects of reinteractions are neglected, the released backward nucleon can 
leave the nucleus keeping its original momentum. These correlated pairs have 
recently been observed in the reaction e + 16 O — > e pp 14 C at low values of 
the energy transfer (180 4 2AQMeV) [45]. 
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To study these correlations it is customary to use the variable a defined as: 



a = (E - P t )/M (5) 

where E and Pi are respectively the energy and longitudinal momentum of 
the Bp and M is the nucleon mass. For Bp, a > 1 since Pi is negative. In 
this model, due to the target motion a correlation between the Bjorken scaling 
variable x and the variable a is expected. In particular the average x for events 
with a Bp is expected to be smaller than in events where no Bp is present, as 
indeed observed in our data (see Table 3). The same correlation is expected 
to hold if the variable v = xy is used, where y = EbadI 'E vis . This variable is 
related to the muon kinematics and can be written as : 

v = (E,- Pl)/M (6) 



where E^ and P 1 ^ are the muon energy and longitudinal momentum. Since, 
generally, in neutrino experiments the muon kinematic variables are well mea- 
sured, v is better suited than x for this purpose. 

The (v , a) correlations were searched for in the data by calculating for each 
event a and v as defined in Eq. 5 and 6. For each a bin we plot the variable 
< vn > defined as: 

< v >Bp 

<v N > = p — (7) 

< V > no Bp 

where < v > no Bp is obtained from the full sample of events without a Bp. 
According to Ref. [10] the average values of v in events with a Bp, < v >Bp, is 
related to the average value of v in events where no Bp is present, < v > no b p , 
by: 

< V > B p = <V > no Bp (2 - a) (8) 



More generally for a cluster composed of £ nucleons the relation is [14]: 

< V > B p = <V> no Bp f 1 - j\ (9) 



Fig. 14a shows < v N > as a function of a. The data indicate a slope of 
—0.22 ± 0.06 with a x 2 / n df — 5.3/6 (the "X 1 jndj in the hypothesis of no de- 
pendence on a is 19.1/7). The two-nucleon correlation mechanism (£ = 2, the 
most probable case when considering the overlapping probabilities of £ nucle- 
ons in the nucleus) fails to describe our data. Either higher order structures 
are playing a leading role [14] or the observed low level of correlation is due to 
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the presence of reinteraction processes. In this case part of the Bp are emitted 
as a result of reinteractions in the nucleus and are not related to the target 
nucleon. The presence of intranuclear cascade processes could therefore dilute 
the existing correlation to the observed level. To test this hypothesis we tried 




a, a 
Fig. 14. The variable < vn > plotted as a function of a. The lines represent the 
predicted correlation (Eq. 9) for a number £ of nucleons in the cluster equal to 2, 
3 and 4- In (a) all the Bp events were used; in (b) only events having a Bp with 
cosOj < 0, 9j being the angle of the Bp with respect to the hadronic jet direction. 

to reduce the component due to rescattering in the selected Bp sample. Hav- 
ing observed the correlation existing between the multiplicities of slow tracks 
and rescattered protons (Sec. 7.2) we applied increasingly tighter cuts on the 
number of slow tracks (P < 700 MeV/ c) . As a consequence of these cuts the 
degree of correlation between a and v increases. The fit values of the slopes 
are reported in Table 6 together with the definitions of the cuts applied. 
We have also observed a strong correlation between the presence of protons 
travelling backward in the lab but forward with respect to the hadronic jet 
direction, and the concentration of events at small Q 2 values and large an- 
gles with respect to the beam. Since also a small Q 2 indicates the presence 
of rescattering, the exclusion of these events should highlight the expected 
correlation. The resulting slope is —0.39 ± 0.07 with a x 2 /ndf = 4.9/6 (see 
Fig. 14b). 

The observed behaviour is consistent with the hypothesis of the correlations 
effects being to some degree hidden by the presence of rescattering. Reducing 
the rescattering component these correlations seem to become stronger. 
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n+ 


n 


slope 


y 2 Indf 

/V I J 


< 3 


< 2 


-0.23 ±0.06 


5.7/6 


< 2 


< 1 


-0.25 ±0.07 


3.3/6 


< 2 





-0.30 ±0.07 


3.7/6 


1 





-0.37 ±0.10 


3.9/6 



Table 6 

The fitted value of the (a,v) slope, and the corresponding \ 2 / n df , for Bp selected 
from events with various numbers of positive (n + ) and negative (n~ ) low momentum 
(P < 700MeV/c) particles. 

8 Conclusions 



We have observed backward proton and 7r _ production in //-nucleus interac- 
tions in the NOMAD detector. 

The slope parameter B of the invariant cross section, parametrized as C e~ Bp2 , 
has been measured and found to be consistent with previous //-nucleus and 
hadron-nucleus experiments. We found that B does not depend on Ehad and 
Q 2 over a wide range of values. This is in agreement with the "nuclear scaling" 
previously observed in hadronic experiments. The observed invariant spectrum 
is not consistent with the existence of a "long" Fermi momentum tail as the 
one proposed in Ref. [33] but agrees with the prediction of Ref. [43]. 
The Bp rate in the NOMAD target (mainly carbon) has been measured and 
compared with the values obtained on different nuclei. While the A dependence 
for neutrino scattering on heavy nuclei is consistent with that of hadron ex- 
periments, and can be parametrized as < Nbp >oc A a with a = 0.68 [21], the 
NOMAD result does not fit this dependence. The A dependence of Bir~ has 
been found to be steeper than that of Bp. 

The backward proton data have been compared with the predictions of reinter- 
actions and short range models. The observed energy dependence is consistent 
with the "formation zone" mechanism. The correlation between the multiplic- 
itiy of slow (P < 700MeV/c) tracks and Bp events indicates the effects of 
reinteractions. However when appropriate cuts are applied in order to reduce 
the intranuclear cascade contributions, the correlation between the Bp and 
the muon scaling variable v, predicted by the short range models, becomes 
stronger. 
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